function y=funacc(x)
y=abs(174.42*(x(1)/x(5))*(x(3)/(x(2)-x(1)))^0.85*sqrt((1-2.62*(1-0.36*(x(4)/x(2))^(-0.56))^1.5*(x(4)/x(2))^1.16)/(x(6)*x(7)))-1.5);